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We generalize the dynamical-mean field theory (DMFT) by including into the DMFT equations depen- 
dence on correlation length of pseudogap fluctuations via additional (momentum dependent) self-energy Ek- 
This self - energy describes non-local dynamical correlations induced by short-ranged collective SDW-like an- 
tiferromagnetic spin (or CDW-like charge) fluctuations. At high enough temperatures these fluctuations can 
be viewed as a quenched Gaussian random field with finite correlation length. This generalized DMFT+Ek 
approach is used for the numerical solution of the weakly doped one-band Hubbard model with repulsive 
Coulomb interaction on a square lattice with nearest and next nearest neighbour hopping. The effective single 
impurity problem is solved by numerical renormalization group (NRG). Both types of strongly correlated 
metals, namely (i) doped Mott insulator and (ii) the case of bandwidth W < U (U — value of local Coulomb 
interaction) are considered. Calculating profiles of spectral densities for different parameters of the model 
we demonstrate the qualitative picture of Fermi surface "destruction" and formation of "Fermi arcs" due to 
pseudogap fluctuations in qualitative agreement with ARPES experiments. "Blurring" of the Fermi surface is 
enhanced with the growth of the Coulomb interaction. 
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Pseudogap formation in the electronic spectrum of 
underdoped copper oxides PE] is is especially striking- 
anomaly of the normal state of high temperature su- 
perconductors. Despite continuing discussions on the 
nature of the pseudogap, we believe that the prefer- 
able "scenario" for its formation is most likely based 
on the model of strong scattering of the charge carriers 
by short-ranged antiferromagnetic (AFM, SDW) spin 
fluctuations E]. In momentum representation this 
scattering transfers momenta of the order of Q = (— , — ) 
(a — lattice constant of two dimensional lattice). This 
leads to the formation of structures in the one-particle 
spectrum, which are precursors of the changes in the 
spectra due to long-range AFM order (period doubling) . 
As a result we obtain non-Fermi liquid like behavior 
(dielectrization) of the spectral density in the vicinity 
of the so called "hot spots" on the Fermi surface, ap- 
pearing at intersections of the Fermi surface with anti- 
ferromagnetic Brillouin zone boundary |z]. 

Within this spin-fluctuation scenario a simplified 
model of the pseudogap state was studied EJ E] un- 
der the assumption that the scattering by dynamic spin 
fluctuations can be reduced for high enough tempera- 
tures to a static Gaussian random field (quenched dis- 
order) of pseudogap fluctuations. These fluctuations are 
characterized by a scattering vector from the vicinity of 
Q, with a width determined by the inverse correlation 
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length of short-range order k = and by appropri- 
ate energy scale A (typically of the order of crossover 
temperature T* to the pseudogap state 0). 

Undoped cuprates are antiferromagnetic Mott insu- 
lators with U 3> W (U — value of local Coulomb in- 
teraction, W — bandwidth of non-interacting band) , so 
that correlation effects are actually very important. It is 
thus clear that the electronic properties of underdoped 
(and probably also optimally doped) cuprates are gov- 
erned by strong electronic correlations too, so that these 
systems are typical strongly correlated metals. Two 
types of correlated metals can be distinguished: (i) the 
doped Mott insulator and (ii) the bandwidth controlled 
correlated metal W ~ U. 

A state of the art tool to describe such correlated 
systems is the dynamical mean-field theory (DMFT) 
[HI 003 HUE]- The characteristic features of correlated 
systems within the DMFT are the formation of inco- 
herent structures, the so-called Hubbard bands, split by 
the Coulomb interaction U, and a quasiparticle (conduc- 
tion) band near the Fermi level dynamically generated 
by the local correlations P 171151 131 ITU] . 

Unfortunately, the DMFT is not useful to the study 
the "antiferromagnetic" scenario of pseudogap forma- 
tion in strongly correlated metals. This is due to the 
basic approximation of the DMFT, which amounts to 
the complete neglect of non-local dynamical correlation 
effects 013 El El EH- As a result, within the standard 
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DMFT approach Fermi surface of a quasiparticle band 
is not renormalized by interactions and just coincides 
with that of the "bare" quasiparticles |7j . Recently we 
have formulated a semiphenomenological DMFT+Ek 
approach JT] , allowing the introduction of a length scale 
(non-local correlations) into DMFT. Below we present 
basic points of this approach with application to the 
Fermi surface renormalization due to pseudogap fluctu- 
ations. 

To include non-local effects, while remaining within 
the usual "impurity analogy" of DMFT, we propose the 
following procedure. To be definite, let us consider a 
standard one-band Hubbard model. The major assump- 
tion of our approach is that the lattice and Matsubara 
"time" Fourier transform of the single-particle Green 
function can be written as: 

Gk(w) = iu; + fi-e{k) - EH -E k H (1) 

where S(ct») is the local contribution to self-energy, sur- 
viving in the DMFT (u = irT(2n + 1)), while E k (a>) 
is some momentum dependent part. We suppose that 
this last contribution is due to either electron interac- 
tions with some "additional" collective modes or order 
parameter fluctuations, or may be due to similar non- 
local contributions within the Hubbard model itself. To 
avoid possible confusion we must stress that £k(*w) can 
also contain local (momentum independent) contribu- 
tion which obviously vanishes in the limit of infinite 
dimensionality d —> oo and is not taken into account 
within the standard DMFT. Due to this fact there is 
no double counting problem within our approach for 
the Hubbard model. It is important to stress that the 
assumed additive form of self-energy £(a>) + Ek(w) im- 
plicitly corresponds to neglect of possible interference of 
these local (DMFT) and non-local contributions. 

The self-consistency equations of our generalized 
DMFT+Ek approach are formulated as follows 

1. Start with some initial guess of local self-energy 
E(w), e.g. E(w) = 0. 

2. Construct Ek(w) within some (approximate) 
scheme, taking into account interactions with 
collective modes or order parameter fluctuations 
which in general can depend on E(w) and /i. 

3. Calculate the local Green function 
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N ^ iuj + n - e(k) - £(w) - E k (w) ' 

(2) 

4. Define the "Weiss field" 

Sq-Vheh + g^h. (3) 



5. Using some "impurity solver" to calculate the 
single-particle Green function for the effective An- 
derson impurity problem, defined by Grassmanian 
integral 

Gd(r-r') =-=?-/ Dc+Dc ta c la (T)c+(T>)eMS cS ) 

^>eS J 

(4) 

with effective action for a fixed ( "impurity" ) i 
r0 r0 

S e $=- / dn / dT2CUr{T\)GQ X {T\ - T 2 )ct (r 2 ) + 

Jo Jo 



+ / drUni^Tjniifr) 
Jo 



Z cff = J DcfDdv exp(-S* cff ), and = T~ l . This 
step produces a new set of values G7 1 (oj). 



6. Define a new local self-energy 

£(w) = Gq 1 ^) - G^{uj). 



(6) 



7. Using this self-energy as "initial" one in step 1, 
continue the procedure until (and if) convergency 
is reached to obtain 



Gu(u>) = G d {u>). 



(J) 



Eventually, we get the desired Green function in the 
form of (I), where E(w) and E k (w) are those appearing 
at the end of our iteration procedure. 

For the momentum dependent part of the single- 
particle self-energy we concentrate on the effects of scat- 
tering of electrons from collective short-range SDW-like 
antiferromagnetic spin (or CDW-like charge) fluctua- 
tions. To calculate Ek(w) for an electron moving in 
the quenched random field of (static) Gaussian spin (or 
charge) fluctuations with dominant scattering momen- 
tum transfers from the vicinity of some characteristic 
vector Q ("hot spots" model [2]), we use the following 
recursion procedure proposed in Refs. |12l which 
takes into account all Feynman diagrams describing the 
scattering of electrons by this random field: 



with 

E n (wk) = A 2 - 



E k (w) = E„ = i(wk) 
s(n) 



(8) 



(5) 



iuj + /i — E(w) — £ n (k) + inv n n — £ n+ i(o;k) 

(9) ' 

The quantity A characterizes the energy scale and 
K = is the inverse correlation length of short 
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range SDW (CDW) fluctuations, e n (k) = e(k + Q) and 
v n = l' L 'k+Q I + l w k+ql f° r n while e„(k) = e(k) and 
v n = \vfc\ + \v^\ for even n. The velocity projections v£ 
and are determined by usual momentum derivatives 
of the "bare" electronic energy dispersion e(k). Finally, 
s(n) represents a combinatorial factor with 



s(n) 



(10) 



for the case of commensurate charge (CDW type) fluc- 
tuations with Q = (7r/a,7r/a) |12| . For incommensurate 
CDW fluctuations J2| one finds 



s(n) 



,J 4r- for odd n 
^ for even n. 



(11) 



If we take into account the (Heisenberg) spin structure 
of interaction with spin fluctuations in "nearly anti- 
ferromagnetic Fermi-liquid" (spin-fermion (SF) model 
Ref. 0j), combinatorics of diagrams becomes more com- 
plicated and factor s(n) acquires the following form 0]: 



s(n) = 



for odd n 
^ for even n. 



(12) 



Obviously, with this procedure we introduce an impor- 
tant length scale £ not present in standard DMFT. 
Physically this scale mimics the effect of short-range 
(SDW or CDW) correlations within fermionic "bath" 
surrounding the effective Anderson impurity. Both pa- 
rameters A and £ can in principle be calculated from 
the microscopic model at hand |11|. 

In the following we will consider both A and espe- 
cially £ as some phenomenological parameters to be de- 
termined from experiments. This makes our approach 
somehow similar in the spirit to Landau approach to 
Fermi-liquids. 

In the following, we discuss a standard one-band 
Hubbard model on a square lattice. With nearest (t) 
and next nearest (t') neighbour hopping integrals the 
"bare" dispersion then reads 

e(k) = — 2i(cos k x a + cos k y a) — At' cos k x a cos k y a , 

(13) 

where a is the lattice constant. The correlations are in- 
troduced by a repulsive local two-particle interaction U. 
We choose as energy scale the nearest neighbour hop- 
ping integral t and as length scale the lattice constant 
a. All energies below are given in units of t. 

For a square lattice the bare bandwidth is W = 8t. 
To study a strongly correlated metallic state obtained 
as doped Mott insulator we use U = AOt as value for the 
Coulomb interaction and a filling n = 0.8 (hole doping). 



The correlated metal in the case of W > U is considered 
for the case of U = At and filling factor n = 0.8 (hole 
doping). For A we choose rather typical values between 
A = O.lt and A = 2i (actually as approximate limit- 
ing values obtained in Ref. 11 ) and for the correlation 
length we have taken £ = 10a (being motivated mainly 
by experimental data for cuprates 0E])- 

The DMFT maps the lattice problem onto an effec- 
tive, self-consistent impurity defined by Eqs. 
In our work we employed as "impurity solver" a reli- 
able method of numerical renormalization group (NRG) 

[num. 

As already discussed in the Introduction, the charac- 
teristic feature of the strongly correlated metallic state 
is the coexistence of lower and upper Hubbard bands 
split by the value of U with a quasiparticle peak at the 
Fermi level. 

Once we get a self-consistent solution of the 
DMFT+Ek equations with non-local fluctuations we 
can compute the spectral functions A(o>,k) for real co: 



1. 



-Im- 
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7r w + fi — e(k) — S(w) — Sk(w) 



> (14) 



where self-energy E(cj) and chemical potential \x are 
calculated self-consistently. Densities of states can be 
calculated integrating (|14|) over the Brillouin zone. 

Extensive calculations of the densities of states, 
spectral densities and ARPES spectra for this model 
were performed in Ref. 11 . In general case pseudogap 
appears in the density of states appears within the in the 
quasiparticle peak (correlated conduction band). Qual- 
itative behaviour of the pseudogap anomalies is similar 
to those for the case of U = 0E], e.g. a decrease of £ 
makes the pseudogap less pronounced, while reducing A 
narrows the pseudogap and also makes it more shallow. 
For the doped Mott-insulator we find that the pseudo- 
gap is remarkably more pronounced for the SDW-like 
fluctuations than for CDW-like fluctuations. Thus be- 
low we present mainly the results obtained using com- 
binatorics (|12(l of spin - fermion model. 

As was noted above within the standard DMFT ap- 
proach Fermi surface is not renormalized by interactions 
and just coincides with that of the "bare" quasiparticles 
jxj- However, in the case of nontrivial momentum de- 
pendence of electron self - energy, important renormal- 
ization of the Fermi surface appears due to pseudogap 
formation 0|. There are a number of ways to define 
Fermi surface in strongly correlated system with pseu- 
dogap fluctuations. In the following we are using in- 
tensity plots (within the Brillouin zone) of the spectral 
density l|14(l taken at u> = 0. These are readily mea- 



sured by ARPES and appropriate peak positions define 
the Fermi surface in the usual Fermi liquid case. 

In Fig. n( a_c ) we show such plots for the case of un- 
correlated metal (U = 0) with pseudogap fluctuations, 
obtained directly from the Green's function defined by 
the recursion procedure ©, ©■ For comparison, in Fig. 
^ (d) we show renormalized Fermi surfaces obtained for 
this model using rather natural definition of the Fermi 
surface as defined by the solution of the equation 

u) - e(k) ReZ(u) - ReE k (w) = (15) 

for to = Q, used e.g. in Ref.jJ- It is seen, that this last 
definition produces the Fermi surfaces close to those de- 
fined by intensity plots of the spectral density only for 
small values of A, while for larger values we can see 
rather unexpected "topological transition" . At the same 
time, spectral density intensity plots clearly demon- 
strate "destruction" of the Fermi surface in the vicinity 
of the "hot spots" with "Fermi arcs" formation with the 
growth of A, similar to those seen in pioneering ARPES 
experiments of Norman et al. |15| . and confirmed later 
in numerous works. 

In Fig. [5] we show our results for the case of corre- 
lated metal with U — it and in Fig. |31 for the doped 
Mott insulator with U = 40i. Again we see the qualita- 
tive behavior clearly demonstrating the "destruction" of 
the well defined Fermi surface in the strongly correlated 
metal with the growth of the pseudogap amplitude A. 
The role of finite U reduces to lower intensity of spectral 
density in comparison with the case of U — and leads 
to additional "blurring" , making "hot spots" less visi- 
ble. Again the "destruction" of the Fermi surface starts 
in the vicinity of "hot spots" for small values of A, but 
almost immediately it disappears in the whole antinodal 
region of the Brillouin zone, while only "Fermi arcs" re- 
main in the nodal region very close to the "bare" Fermi 
surface. These results give a natural explanation of the 
observed behavior and also of the fact that the existence 
of "hot spots" regions was observed only in some rare 
cases [IE] . 

For the case of doped Mott insulator (U — 40t) 
shown in Fig.[5]we see that the "Fermi surface" is rather 
poorly defined for all values of A, as the spectral den- 
sity profiles are much more "blurred" than in the case 
of smaller values of U, reflecting important role of cor- 
relations. 

It is interesting to note that from Figs. [3 El it is 
clearly seen that the "natural" definition of the Fermi 
surface from Eq. I|15[) is quite inadequate for correlated 
systems with finite U and nonlocal interactions (pseudo- 
gap fluctuations) , signifying the increased role of strong 
correlations. 



To summarize, we propose a generalized DMFT+Sk 
approach, which is meant to take into account the im- 
portant effects due to non-local correlations in a system- 
atic, but to some extent phenomenological fashion. The 
main idea of this extension is to stay within a usual effec- 
tive Anderson impurity analogy, and introduce length 
scale dependence due to non-local correlation via the 
effective medium ("bath") appearing in the standard 
DMFT. This becomes possible by incorporating scatter- 
ing processes of fermions in the "bath" from non-local 
collective SDW-like antiferromagnetic spin (or CDW- 
like charge) fluctuations. Such a generalization of the 
DMFT allows one to overcome the well-known short- 
coming of k-independence of self-energy of the stan- 
dard DMFT. It in turn opens the possibility to access 
the physics of low-dimensional strongly correlated sys- 
tems, where different types of spatial fluctuations (e.g. 
of some order parameter), become important in a non- 
perturbative way at least with respect to the important 
local dynamical correlations. However, we must stress 
that our procedure in no way introduces any kind of sys- 
tematic 1/d-expansion, being only a qualitative method 
to include length scale into DMFT. 

In our present study we addressed the problem of the 
Fermi surface renormalization ("destruction") by pseu- 
dogap fluctuations in the strongly correlated metallic 
state. Our generalization of DMFT leads to non-trivial 
and in our opinion physically sensible k-dependence of 
spectral functions, leading to Fermi surface renormal- 
ization quite similar to that observed in ARPES exper- 
iments. 

Similar results were obtained in recent years using 
the cluster mean-field theories ^7j. The major advan- 
tage of our approach over these cluster mean-field the- 
ories is, that we stay in an effective single-impurity pic- 
ture. This means that our approach is computationally 
much less expensive and therefore easily generalizable 
for the account of additional interactions. 
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Fig.l. Fermi surfaces obtained for uncorrelated case of 
[7 = and filling factor n = 0.8. Shown are intensity 
plots of spectral density 11411 for ui = 0. (a) - A = O.lt; 
(b) - A = 0.3t; (c) - A = t; (d) - "Fermi surfaces" 
obtained solving Eq. 1151 . Dashed line denotes "bare" 
Fermi surface. 




Fig. 2. "Destruction" of the Fermi surface as obtained 
from the DMFT+E k calculations for U = 4i and 
n = 0.8. Notations are the same as used in Fig. 
(a) - A = 0.2t; (b) - A = OAt; (c) - A = t; (d) - 
A = 2t. Black lines show the solution of Eq. 1151 . 
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